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We study the folding of the regular triangular lattice in three dimensional embedding 
space, a model for the crumpling of polymerised membranes. We consider a discrete 
model, where folds are either planar or form the angles of a regular octahedron. These 
"octahedral" folding rules correspond simply to a discretisation of the 3d embedding space 
as a Face Centred Cubic lattice. The model is shown to be equivalent to a 96-vertex model 
on the triangular lattice. The folding entropy per triangle In q^d is evaluated numerically 
to be q^d = 1.43(1). Various exact bounds on q 3 d are derived. 
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Polymerised membranes have recently been the object of intense investigation, both 
as natural 2d generalisations of polymers and as models for biological systems. From the 
theoretical point of view, polymerised membrane models are archetypical examples of the 
interplay between geometry and statistical mechanics. In this paper, we consider only 
phantom membranes, where steric constraints due to self-avoidance are not taken into ac- 
count. As in other domains of random surface physics, the recourse to discrete formulations 
provides a powerful approach, with the advantage of permitting direct numerical simula- 
tions, as well as leading to exact solutions for particular models. In the simplest discretised 
version, the polymerised membrane is modelled by a regular triangular lattice with fixed 
connectivity, embedded in a d-dimensional space. In one class of models, the lengths of the 
links of the lattice are allowed to have small variations. These tethered membranes were 
first studied in |1J , where a geometrical 3d continuous crumpling transition was predicted. 
This prediction was corroborated by various analytic results on continuous models 
and other numerical simulations ||. In a second class of models, the lengths of the links 
of the lattice are fixed, say to unity. For such membranes, the only remaining degree of 
freedom is that of folding of the lattice, whose links serve as hinges between neighbouring 
triangles. Membrane folding was first studied in 2d embedding space, where it can be 
formulated as an 11-vertex model on the triangular lattice. The entropy of folding, which 
counts the number of distinct folded states of the lattice in the thermodynamic limit, was 
first estimated numerically in ||. The 2d folding problem was then solved exactly in 0, 
through the equivalence with the 3-colouring problem of the links of the triangular lat- 
tice ||. A further study, including a bending rigidity, performed in Pj, led to numerical 
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evidence of a first order folding transition of the membrane. 

In this paper, we address the 3d folding problem, where the embedding space now has 
3 dimensions, in an attempt to recover results for tethered membranes in the folding lan- 
guage. In particular the 2d folding transition could be smoothed and lead to a continuous 
3d folding transition, analogous to the crumpling transition. The present work is devoted 
to a suitable definition of the 3d folding problem, amenable to numerical simulations, and 
to various reformulations as a vertex model, constrained spin system and colouring prob- 
lem. We use these formulations to estimate the entropy of 3d folding both numerically and 
analytically. 

The paper is organised as follows. In the first section, we define the general problem 
of folding of the triangular lattice in embedding space with arbitrary dimension and recall 
a few results for the case of dimension 2. In section 2, we define the discrete problem 
of 3d "octahedral" folding, which corresponds to a discretisation of the 3d embedding 
space as a Face Centred Cubic lattice. In section 3, we show that the 3d octahedral 
folding is equivalent to a 96-vertex model on the triangular lattice. This is done by use 
of face spin variables defined in section (3.1), which in turn parametrise some generators 
of the tetrahedral group studied in section (3.2). From the group formulation we derive 
in section (3.3) the local folding constraints on these face variables, leading to a complete 
determination of the 96 folding vertices. Section 4 is devoted to the numerical evaluation 
of the entropy of 3d folding. Various exact bounds on this entropy are derived in section 
5, by reexpressing the 3d octahedral folding problem as a dressed 3-co louring problem 
(section (5.1)), and performing some rough bounding of the partition function (section 
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(5.2)). Improved bounds are obtained in section (5.3) by comparison with the 2d folding 
problem in a field. We discuss in section 6 our results for the entropy (section (6.1)), as well 
as a natural generalisation of the octahedral folding problem in dimension d > 3 (section 
(6.2)). Concluding remarks are gathered in section 7, and some extensions of section (5.2) 
are given in the Appendix. 

1. The folding rule 

A folding in R d of the regular triangular lattice is a mapping which assigns to each 
vertex v of the triangular lattice a position X„ in <i-dimensional embedding space R d , with 
the "metric" constraint that the Euclidean distance |X„ 2 — X Wl | in R d between nearest 
neighbours v\ and V2 on the lattice is always unity. Under such a mapping, each elementary 
triangle of the lattice is mapped onto an equilateral triangle in R d . In general, two adjacent 
triangles form some angle in R d , i.e. links serve as hinges between triangles and may be 
(partially) folded. 




Fig. 1: The oriented triangular lattice: triangles pointing up (resp. down) 
are oriented counterclockwise (resp. clockwise). The three tangent vectors 
tj, i = 1,2, 3, have a vanishing sum in the embedding space. 



Folding is best described in terms of tangent vectors, which are link variables defined as 
follows: we first orient the links of the lattice as in Fig.l, with triangles pointing up (resp. 
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down) oriented counterclockwise (resp. clockwise), and define the tangent vector between 
two neighbours V\ and i>2 as the vector 

if the arrow points from v± to i>2. The metric constraint states that all tangent vectors 
have unit length. Moreover, with our choice of orientation, the three tangent vectors tj, 
i = 1,2,3, around each face of the lattice must have vanishing sum. This is the basic 
foldins; ru lJ: 

ti+t 2 + t 3 = 0. (1.2) 

Up to a global translation in R d , a folding is therefore a configuration of unit tangent 
vectors defined on the links of the lattice, obeying the folding rule ( |1 . 2j ) around each 
triangle. 

The two-dimensional (d = 2) folding problem of the triangular lattice was addressed 
in 0. It is easy to check that, up to a global rotation in the embedding plane, all the 
link variables are forced to take their value among a fixed set of three unit vectors with 
vanishing sum. This permits a reformulation of the d = 2 folding problem as that of the 
3-colouring of the links of the triangular lattice: calling the three fixed vectors blue, white 
and red, the folding rule translates into the constraint that the three colours around each 
triangle have to be distinct. This 3-colouring problem was solved by Baxter in || by use 
of Bethe Ansatz techniques. His result for the thermodynamic partition function measures 

1 On the dual hexagonal lattice, whose vertices are at the centers of the triangles, this basic 
rule translates into a local conservation law at each vertex, reminiscent of the so-called "ice rule" 
of integrable vertex models. 
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the number of folding configurations Z^d oc q 2d A for a lattice with triangles, in the 
limit of large This gives the folding entropy per triangle s 2 d = logfed), with 

q 2d = ^r(l/3) 3 / 2 = 1.20872... (1.3) 

The 2d folding problem has also been studied in the presence of bending rigidity, which 
associates an energy to each folded link, and a magnetic field coupled to the normal vector 
to the triangles 0. The system was found to undergo a first order folding transition. At 
zero rigidity and zero magnetic field, the lattice is in an entropic folded phase. At large 
enough rigidity and/or magnetic field, the lattice becomes totally unfolded. 

In this paper, we address the case of embedding in R d=3 and, in particular, we are 
interested in computing the folding entropy s^d- 



2. Octahedral Folding 



In the general 3-dimensional folding problem, the local folding constraint (|1.2|) imposes 
only that the three tangent vectors around each face be in the same plane and have relative 
angles of 27r/3. This, however, does not impose any constraint on the relative positions 
of the two planes corresponding to two adjacent faces, which may form some arbitrary 
continuous angle. As opposed to the 2d case, this then leads to a problem with continuous 
degrees of freedom. 

Here we wish instead to define some discrete model of folding in 3<i, in which only a 
finite number of relative angles are allowed between adjacent faces. More specifically, we 
shall also impose that the link variables themselves take their values among a finite set of 
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Fig. 2: The oriented octahedron: the edges around each face form triplets of 
tangent vectors with vanishing sum. The four normal vectors n^, i = 1, 2, 3, 4, 
are represented on the corresponding outward oriented faces. 

tangent vectors, now in R 3 . For symmetry reasons, we will take this set of tangent vectors 
to be the (oriented) edges of a regular solid of R 3 , made of equilateral triangles only. 

There are only three regular solids in R 3 made of equilateral triangles: the tetrahedron, 
the octahedron and the icosahedron. The edges of the tetrahedron (resp. icosahedron) , 
however, cannot be consistently oriented in order for the corresponding tangent vectors 
to satisfy ( fT2|) around each face. This is because each vertex is surrounded by an odd 
number 3 (resp. 5) of triangles. There is no such problem for the octahedron, as shown 
in Fig. 2. The 12 links of the octahedron are oriented consistently to form 8 triplets of 
tangent vectors with vanishing sum corresponding to the 8 faces of the octahedron. From 
now on, we shall therefore consider the restricted 3d "octahedral folding" problem, where 
the tangent vectors are chosen from the set of the 12 edge vectors of a regular oriented 
octahedron. In the folding process, the folding rule ( |1.2| ) imposes that the three links of a 
given face on the original triangular lattice are mapped onto one of the 8 triplets of tangent 
vectors above. For a given triplet, the triangle can still be in 3! states corresponding to 
the 3! permutations of the three edges. Each triangle can therefore be in one of 48 = 8 x 6 
states. 
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The 8 faces of the octahedron can be labelled as follows: we consider for each face 
its normal vector, pointing outwards or inwards according to the orientation of its tangent 
vectors on the octahedron (see Fig. 2). There are four outward oriented and four inward 
oriented faces which alternate on the octahedron. The normal vectors to opposite faces 
are equal. This thus defines a set of only four vectors ni, 112, 113 and 114, furthermore 
satisfying the sum rule ri! + n 2 + n 3 + n 4 = 0. Each face is labelled by its orientation 
(outward or inward) and its normal vector (1, 2, 3 or 4). Notice that, by labelling the 
normal vectors, we can arbitrarily fix one of the two possible chiralities of the octahedron 
in this construction. We choose the one corresponding to Fig. 2. 



Fig. 3: The Face Centred Cubic lattice viewed as a packing of 3d space with 
octahedra and tetrahedra. 

The 12 oriented edge vectors of the octahedron are actually identical to the 6 edge 
vectors of a tetrahedron, now taken with the two possible orientations. The four normal 
vectors above are also the normals to this tetrahedron. For each folding map, the image 
of the folded lattice in R 3 lies therefore on a 3d Face Centred Cubic (FCC) lattice, which 
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consists of a filling of space by octahedra complemented by tetrahedra, as shown in Fig. 3. 
In this respect, the "octahedral folding" problem simply corresponds to discretising the 
embedding space as a FCC lattice. 

Note finally that nothing prevents the lattice from intersecting itself, hence our con- 
struction describes a phantom membrane. The introduction of self-avoidance would result 
in much more elaborate non-local constraints, far beyond the scope of the present study. 

3. 96-vertex Model 

When stated in terms of tangent vectors, the 3d "octahedral folding" problem involves 
three types of constraints: face, link and vertex constraints. The first constraint, around 
each face, imposes that the three tangent vectors of a given triangle form one of the 8x6 
(ordered) triplets with vanishing sum. The second constraint, on each link, arises because 
two adjacent triangles share a common tangent vector. Given the state of one triangle, 
any adjacent triangle has one of its tangent vectors already fixed and thus is left with only 
4 = 48/12 possible states. 



acute 




Fig. 4: The four possible folding angles between two adjacent triangles. The 
neighbour of the dark triangle may (i) be itself on top of the dark triangle 
(complete fold), (ii) occupy the symmetric position in the same plane (no 
fold), (iii) lie on the same octahedron (i.e. form an obtuse angle) or (iv) lie 
on the same tetrahedron (i.e. form an acute angle). 
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They correspond simply to the four values for the relative angle between two neighbouring 
triangles, i.e. the angle between the normal vectors, depicted in Fig. 4. These four values 
are (no fold: the triangles are side by side), 180° (complete fold: the triangles are on 
top of each other), arccos(— 1/3) ~ 109°28' (fold with acute angle: the two triangles lie on 
the same tetrahedron) and arccos(l/3) ~ 70°32' (fold with obtuse angle: the triangles lie 
on the same octahedron). Finally, there is a third constraint on the six successive folds 
around each vertex of the lattice: after making one loop, the same tangent vector must be 
recovered. Since the "metric constraint" is local, there are actually no constraints other 
than these three (face, link and vertex) constraints. 

In the study of the 2d folding problem, i.e. of the 3-colouring problem, the face and 
link constraints are taken into account by going to Z2 spin variables Oi defined on the 
faces of the lattice. Ordering the colours cyclically, the spin is +1 (resp. —1) if the colour 
increases (resp. decreases) from one link to the neighbouring one on the triangle, oriented 
counterclockwise. In this language, the actual folds take place exactly on the domain walls 
of the spin variable. Instead of having a Z3 colour variable per link, one is left with a Z2 
spin variable per triangle. The vertex constraint translates into a constraint on the six 
spins o"i, ...,06 around each vertex of the lattice, namely that Y^i=i °"i = mod 3. This 
leads to 22 possible local spin configurations around each vertex, or equivalently, after 
removing the global Z2 degeneracy of reversal of all spins, to an 11-vertex model on the 
lattice 0. 

In this section, we shall proceed in the same way for the 3d "octahedral folding" and 
account for the face and link constraints by expressing folded configurations in terms of two 
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Z 2 variables on the triangles. These variables will indicate the relative states of successive 
links around the face. We shall then count the number of allowed hexagonal configurations 
around a vertex: we will find that our problem is now expressible as a 96-vertex model. 
These vertices and the corresponding rules on the Z 2 variables will be identified in the 
next section. 

3.1. Face variables and counting of the vertices 

© n i 

► Co) 




Fig. 5: The labelling of an edge of the octahedron, 1 < % 7^ j < 4, 
according to the adjacent normal vectors. 

Let us label the 12 edges of the octahedron as follows: each edge is shared by two 

adjacent faces, one outward and one inward oriented (see Fig. 2). We label the edges by 

the indices 1 < i 7^ j < 4, when the normal vector for the outward face is and the 

one for the inward face is nj. This convention is illustrated in Fig. 5. There are 12 such 

couples The unit tangent vector associated with the link (z, j) is given by 

3 

= ^ n * x (~ n j)- (3- 1 ) 

Consider now an elementary triangle of the lattice. Starting from one of its links 

the subsequent link (k, I) counterclockwise must share a face with on the octahedron. 

This leads to the 4 following possibilities, labelled by the Z 2 face variables z, a E {±1}: 
z = +l: (i,j) -> (i,Z), / ^ j, tiji = -a = ±1 

(3.2) 

z = -l: -> (k,j), k^i, e ijk = +a = ±1 
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Fig. 6: The transition from a link (i, j) to a subsequent link (k, I) is described 
by the two Z 2 face variables z and a. 

where = J2i e ijki is defined in terms of the totally antisymmetric tensor e^, equal 

to the signature of the permutation (ijkl) of (1234). The value z = +1 (resp. z = —1) 

indicates that the two tangent vectors share an outward oriented (resp. inward oriented) 

face on the octahedron. The spin variable a takes the value +1 (resp. —1) if (fe, I) follows 



(resp. precedes) on their common (oriented) face of the octahedron. Using (|3.1| ), one 
can check that the variable a also indicates whether the normal vector to the triangle (in 
the embedding space R 3 ) is parallel (cr = +1) or antiparallel (a = — 1) to the corresponding 
normal vector of the octahedron. 

Considering now two neighbouring triangles, the 4 possible relative values z%/zi and 
cr 2 /c7i indicate which type of fold they form, with the correspondence displayed in Table I 
below. 



Z2/Z1 


02M 


angle 


1 


1 


no fold 


1 


-1 


complete fold 


-1 


1 


acute fold 


-1 


-1 


obtuse fold 



Table I: The relative folding state of two neighbouring triangles according 
to their relative values of z and a. 

The domain walls for the z variable are the location of the folds which are either acute 

or obtuse, whereas those for the a variable are the location of the folds which are either 
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complete or obtuse. The superposition of these two types of domain walls fixes the folding 
state of all the links, specifying the folding state of the lattice up to a global orientation. 

The use of z and a variables instead of the 12 variables incorporates the face and 
link constraints. Like in the 2d case, the vertex constraint is more subtle and will be studied 
in the next section. Nevertheless, we can easily count at this stage the number of possible 
configurations around a vertex satisfying this constraint, i.e. the number of possible folded 
states of an elementary hexagon. Indeed the mapping ( |3.2| ) may be represented by a 12 x 12 
connectivity matrix Mujy^a with i ^ j and k ^ I: 

M(i,j),(k,i) = &ik + Sji — 2 5 ik 5ji . (3.3) 

This matrix acts as a transfer matrix between two successive internal links of the hexagon. 
The number of configurations of a hexagon is simply given by: 

Tr(M 6 ) = 4608 , (3.4) 

where the trace guarantees that the same link variable is recovered after one loop. These 
4608 configurations count as distinct all the foldings which are related by a global change 
of orientation of the hexagon in embedding space. The order of the resultant degeneracy 
is 48, corresponding to 12 choices for the first tangent on the octahedron times 4 for the 
choice of the second from among its 4 neighbours (this latter choice corresponds to the 
4 choices of the z and a variables on the corresponding triangle). This leaves us with 
4608/48 = 96 distinct configurations. 

The above computation is also equivalent to counting the number of closed paths of 
length 6 on the cuboctahedron, i.e. the solid whose vertices sit at the centers of the edges 
of the octahedron, and whose connectivity matrix is given by ( |3.3j ). 
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3.2. The tetrahedral group 

In order to derive the vertex constraints for the z and a variables around any vertex 

of the lattice, we will use a group formulation of (|3.2|) . We note here that the 12 links 

of the octahedron may be uniquely labelled by the even permutations of the set 1, 2, 3, 4. 

Indeed the link is equivalently represented by the even permutation (ijkl), with 

tijki = 1- With this labelling the alternating group A4, the group of even permutations of 

four elements, also known as the tetrahedral group, acts on the space of links. In other 

words, each element of is a one-to-one mapping of the set of links onto itself. A fixed 

element of A4 acts on a given link by permutation of the four link labels. The order of A4 

is 4!/2 = 12. These 12 elements map a given link to exactly the 12 links of the octahedron. 

The group A4 can be generated by two elements 7*1 and T3 defined as: 

n = (243) 

(3.5) 

r 3 = (134) 




Fig. 7 : The action of the generators n and T3 of the tetrahedral group A4 on 
a particular edge of the oriented octahedron. All the edges of the octahedron 
are equivalent, and n and T3 act simultaneously on all the edges according 
to the same picture. 

and corresponding to the two elementary mappings shown in Fig. 7. Here we use the 
standard notation for permutations, by listing their cycles not reduced to one element 
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(e.g. (243) is the permutation (1234) — > (1423), whereas (12) (34) is the permutation 
(1234) -> (2143)). 

In terms of these generators the 12 elements of A4 are: 



e 


3 

= T{ = 


3 • i j ■ j 

r 3 = identity 




= (243) 




T 3 


= (134) 




2 

T2 = r x 


= (234) 




2 

r 4 = r 3 


= (143) 




T1T3 


= (132) 




T3T1 


= (124) 




^22 

r l r 3 


= (142) 




^2 2 
r 3 r l 


= (123) 




9 = r i r 3 2 


= 7"3 r\ - 


= (14) (23) 


d = r^r 3 


= r h\ -- 


= (13) (24) 


/ = T1T3T1 


= T3T1T3 


= (12)(34). 



(3.6) 



Any other sequence of t\ and r 3 can be reduced to one of the above elements by use of the 
four relations! 

7*! = r 3 = identity 
rfra = rfn (3.7) 
T1T3T1 = T3T1T3 • 

These relations can be easily understood graphically by following the successive images of 

a link around the octahedron (see Fig. 7 above). 
2 These relations also imply that rir| = T3T1 . 
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3.3. Vertex rules 

Starting from the link, say the four choices of subsequent link (|3.2| ) correspond 

to the application of the four operators ti, T2, T3 or T4, with the correspondence given in 
Table II below. 



z 


a 


group element 


1 


1 


T\ 


1 


-1 


T 2 


-1 


1 


T3 


-1 


-1 


T4 



Table II: The correspondence between the z and a variables and the elements 
r of A 4 . 

A folding of an elementary hexagon of the regular triangular lattice corresponds to a 
product of six basic group elements r a , a G {1,2,3,4}, chosen from the four operators 
above, and such that 




Fig. 8: The six Z{ and Oi variables around a given vertex, and the colours Cj 
of the interior links. 

Let us now translate this constraint into folding rules on the six z and a variables 
around each vertex (see Fig. 8). A first folding rule involves the a variable alone. It 
ensures the 3-colourability of the links of the triangular lattice in the following way. Let 
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us assign one of three colours (0=blue (B), l=white (W) or 2=red (R)) to each link of 
the octahedron. Choosing the colour of link (1342) to be say, (blue), the colour of a link 
(ijkl) is obtained by counting the total number of t\ and T3 operators required to reach 
this link from (1342). The colour is this number modulo 3. Since the use of relations (|3.7|) 
preserves the total number of t\ plus T3 mod 3, the colour is well-defined. 




Fig. 9: The assignment of colours for the edges of the octahedron. The four 
edges of a given colour lie in the same plane. The three colours around a face 
are distinct. 



This is illustrated in Fig. 9: there are 4 blue (resp. white, red) edges lying in the same 
plane; the three edges have distinct colour around each face. Each folding of the triangular 
lattice induces therefore a 3-colouring of its links. Conversely, a given 3-colouring of the 
triangular lattice does not specify its 3d folding state entirely since the same colour may 
correspond to 4 distinct edges of the octahedron. 

The 3-colourability requirement still leads to a first constraint on the a variables 
around each vertex: the total number of n and T3 mod 3 in the product P (where we 
substitute ri = rf an d T4 = t|) must be a multiple of 3. Now in the assignment of colours 
7~i and T3 count for 1, while T2 = rf and T4 = r| count for 2 = — 1 mod 3. Since this 
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corresponds precisely to their respective values of a, we find the first folding rule: 
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^Vi = 0mod3 (3.9) 



for the six spins around the central vertex of the hexagon. As mentioned before, this rule 
first emerged in the 2d folding problem in J7| , where it also guaranteed the 3-colourability 
of the triangular lattice. 

In contrast with the 2d situation, the restriction ( |3.9|) is not the only constraint here. 
It only ensures that the product P is one of the four colour-preserving elements e, d, g or 



/ of (|3lf). On the other hand, any sequence of t\ and Ts satisfying (|3.9| ) may be naturally 

written as a product of e, d, g and / operators by simply regrouping the r's into triplets. 

It can be checked that e, d, g and / form a Z2 x Z2 subgroup of A4 with the representation 

e = rf = rf = (1,1) 

g = riT | = t z tI = (77,1) 

(3.10) 

d = rfr 3 = rfri = (1, x) 

f = tit 3 ti = T3T1T3 = (77, x), 
where rf 1 = \ 2 = 1. For the product P to be the identity e, we need the number of r\ 

and the number of x i n this decomposition to be separately even. Consider the product P 

written as a sequence of n and T3 operators and mark the spaces between r's by 0, 1 or 2 

consecutively: 

210210210210 

lillillJIIII 

TTTTTTTTTTTT 



The number of 77 (resp. x) is the number of changes from t\ to T3 or T3 to n occurring 
at the spaces labelled by 2 (resp. 1). This follows directly from ( |3.10| ). On the hexagon 
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these changes can occur only at the border between two neighbouring triangles i.e. on 
interior links. The position of a change (0, 1 or 2) is then simply the colour of the link at 
which it occurs. In terms of a variables, the colour c\ of the internal link i 7 i = 1, 6, is 
given by 

i 

Ci = co + ^^Oj mod 3 (3-H) 
i=i 

(the first folding rule ( |3.9| ) ensures that cq = cq). The quantity 1 ~ Zl 2 2l+1 is or 1 depending 
on whether or not a change occurs at link i (with the convention Zf = Z\). Defining 

6 i_ z . . 

a ° = p^(5(c„cmod 3) (3.12) 

i=l 

as the total number of changes on the links of colour c, a c must be even. This gives the 
second folding rule 

a c = mod 2 : c = 1,2 . (3.13) 

Note that ( |3.13| ) implies that «o = mod 2 since cto + ct\ + «2 = mod 2, as the total 
number of changes between t\ and is even by cyclicity. 

To understand this second folding rule, we first note that a change from n to T3, or 
vice versa, corresponds to a switch to another face of the octahedron, therefore crossing one 
of the blue, white or red edge planes bisecting the octahedron (see Fig. 9). The requirement 
of returning to the same (say blue) link after 6 steps is then equivalent to that of crossing 
each of the white and red planes an even number of times. Indeed a.\ (resp. 0,2) is simply 
the number of crossings of the white (resp. red) plane. 

The two folding rules ( jOl) and ( |3.13|) are equivalent to the condition (|3j8|) and there- 
fore characterise the vertex constraint entirely. 
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Fig. 10: The 96 vertices satisfying the two folding rules (|3.9[ ) and ( |3.13| ): 
no line corresponds to no fold, a thick line corresponds to a complete fold, 
a thin line corresponds to a fold with obtuse angle and a dashed line corre- 
sponds to a fold with acute angle. The degeneracy of each vertex under cyclic 
permutations of the links is indicated. 

With the two folding rules ( |3.9j ) and ( |3.13j ) we find 384 = 96 x 4 vertex configurations 
(there is a 4-fold global degeneracy under reversal of z or u), as predicted by our previous 
counting. The 96 folding vertices are displayed in Fig. 10 with the following conventions: 
no line corresponds to no fold; a thick line corresponds to a complete folding (180°, flip of 
a only); a thin line corresponds to a fold with obtuse angle (arccos(l/3) ~ 70°32' between 
normal vectors, flip of both a and z) and finally a dashed line corresponds to a fold with 
acute angle (arccos(— 1/3) ~ 109°28', flip of z only). The degeneracy of each vertex under 
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Fig. 11: Examples of 3d octahedral foldings of an elementary hexagon and 
the corresponding vertices of Fig. 10. 

cyclic permutations of the links is also indicated. 

In Fig. 11 we display a few examples of vertices and the corresponding foldings in 3d space. 

The 3d octahedral folding problem is now represented as a 96-vertex model with the 
vertices of Fig. 10. 

4. Entropy: numerical calculation 

In this section, we use the above folding rules to estimate numerically the 3d folding 
entropy per triangle s 3d = In <?3d- 

The partition function for a parallelepiped of 2L x M triangles, also conveniently 
thought of as a rectangle ofLxM squares divided into two triangles along one diagonal, 
can be written as 

Z LM = Tr T(L) M (4.1) 

in terms of a row to row transfer matrix T(L). 

The above reexpression of the folding problem in terms of two face variables, the spins 
z, o = ±1, enables us to write the transfer matrix from a row of 2L triangles to another 
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as a 16 L x 16 L matrix (there are 4 spin configurations (z, a) on each of the 2L triangles). 
We choose for instance to apply free boundary conditions on the right and left sides of the 
rectangle. The entries of the row to row transfer matrix read 



, z '° / 

°0 / Z'l 

/ °' 


°2 / Z' 3 






, z 2L-2/ 
°2L-2 / 

X Z 2L-1 
/ C 2L-1 


z / 

° / z 
/ 1 


*- i* 

z 2 / 

/ Z 3 
/ CT 3 


' 1,2* 


L-l 1 


z 2L-2/ 
°2L-2 / 

/ z 2L-1 
/°2L-1 



(4.2) 



n u i x 



where Uj and V} respectively impose the first and second folding rules ( ]3.9| )-( |3TT5| ) around 
each inner vertex j, namely 

t/j = 5(a 2 j-2 + &2j-l + &2j + &2j-l + a 2j + a 2j+li® m °d 3) 
2 

^' = II S ( a c( Z 2ji Z 2j-l: z 2j-2, z 2j-l, 4j: 4j+l)> mod 2 ) 



(4.3) 



c=l 



with Q! c (-Zi, -Z2, -23, £4, Z5, ^e) defined by eq.( |3.12|) . Note that Vj also depends implicitly on 
the cr's around the vertex j. 

In the large M limit, the partition function ( |4 . 1| ) is dominated by the largest eigenvalue 
Xl of the transfer matrix T(L) and the corresponding free energy per row reads 



lim — In Zt m 



In Aj 



(4.4) 



The thermodynamic entropy per triangle of the 3d folding problem is therefore 

s 3d = lim —-pin A = In q 3d . 

L — >oo ZIj 



1 



(4.5) 
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The numerical calculation of A^ is simplified by the fact that T(L) is a sparse matrix, 
with entries or 1. Actually the total number of non-vanishing elements of T(L) reads 

N L = T ( L ){^},{z',*'} = Tr(j[T(2)} L ) (4.6) 

{<r,z},{a',z>} 

where J is the 256 x 256 matrix with all entries equal to 1 (see ||). For large L, Nl 

is dominated by the Lth power of the largest eigenvalue A2 of T(2), easily computed 
numerically (see Table III below), henceS 

N L oc (25.48912. ..) L (4.7) 

to be compared with the total number of elements of T(L), (256) L . 

We display the results for the number of non-vanishing elements and the largest 
eigenvalue A^ in Table III below. 



L 


N L 


A L 


(A L /A L _!)5 


1 


256 


16.0000000000 




2 


6144 


25.4891252930 


1.2621688994 


3 


153600 


44.9037351935 


1.3272837214 


4 


3891200 


83.9628811670 


1.3674215892 


5 


98992128 


162.6428871867 


1.3917905107 


6 


2521694208 


322.3383026700 


1.4077917628 



Table III: The number iV^ of non-vanishing transfer matrix elements and 
the maximum eigenvalue Xl for strips of width L = 1,2,3, ...,6. The ratio 
(Al/Al-i)^ is less sensitive to finite size effects than A^ 2L . 

The entries of the matrix T(L) are generated by gluing together around a central 

vertex two transfer matrices of respective sizes L/2,L/2 if L is even or (L — l)/2, (L + 1)/2 

3 The largest eigenvalue of T(2) matches the value 14 + y/132 up to the precision of our 
calculation. 
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if L is odd. In this way, the time used for generating T(L) becomes comparable to that 
used for extracting its largest eigenvalue, by iterative applications of T(L) on a vector, 
which we normalise at each step (the normalisation factor converges to The value 

L = 6 is obtained by applying some extra symmetry arguments to reduce the size of the 
matrix T(6). 

Using the Pade-Shanks transformation |T^,JTJ, we extrapolate the values of (A^/A^-i) 1 / 2 



of Table III above to get a numerical estimate of the partition function per triangle Q4.5Q . 
We find 

q 3d ~ 1.43(1) . (4.8) 



5. Dressed 3— colouring and various bounds on the entropy 

In this section, we derive exact bounds for the partition function per triangle qsd- For 
this purpose, we turn back to the initial definition in terms of tangent vector link variables. 

5.1. 3d folding as dressed 3-colouring 

The link variable of the 3d folding problem takes its values among the 12 tangent 
vectors to the octahedron of Fig. 2. Let ei, e2, e3 denote the canonical basis of R 3 and let 
us fix the positions of the 6 vertices of the octahedron to be (±l/y/2, 0, 0), (0, ±l/v2, 0) 
and (0, 0, ±l/y2) in this basis. For instance, the face with normal vector ni pointing out 
in Fig.2 has the vertices (1^/2,0,0), (0,1/^,0), (0,0,1/^). Denoting by x, = e,/V2, 
i = 1,2, 3, we list in Table IV below the 12 unit tangent vectors grouped according to their 
outward oriented face label. We also display the colour of each vector, as defined in Fig. 9. 
In the present language, the colour of a tangent vector simply corresponds to the index of 
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the missing basis vector x modulo 3, i.e. 3 = 0^ B, 1 — > W and 2 — > i?. 



face 


normal vector 


tangent 


colour 






-Xi + x 2 


B 


1 


ni = ^(xi + x 2 + x 3 ) 


-x 2 + x 3 


W 






-x 3 + Xi 


R 






-xi - x 2 


B 


2 


n 2 = y^|(xi - x 2 - x 3 ) 


x 2 - x 3 


W 






x 3 + Xi 


R 






xi - x 2 


B 


3 


n 3 = \fi(-*i ~ x 2 +x 3 ) 


x 2 + x 3 


W 






-x 3 - Xi 


R 






Xl + x 2 


B 


4 


n 4 = ^(-xi + x 2 - x 3 ) 


-x 2 - x 3 


W 






x 3 - Xi 


R 



Table IV: The twelve unit tangent vectors around the faces of the octahedron 
of Fig. 2. We display the face label, the normal vector to the face (pointing 
outwards), the three tangent vectors around the face and their respective 
colour. 



8 iXi + B jXj 




Fig. 12: A typical configuration of tangent vectors around a triangle, ex- 
pressed in the x basis. Once the colour of each link is chosen, namely the 
three couples (xj,Xj), (xj,Xfc) and (x/-,Xj), the three signs e^, ej and efc are 
still arbitrary. 



The 48 configurations of tangent vectors around a face of the triangular lattice are 
specified by the following two sets of data, displayed in Fig.12: 
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(i) one of the 6 allowed colourings of the three links, corresponding to a permutation 
of B,W, R. This specifies for each link the plane Xj, Xj where the tangent vector lies. 

(ii) the assignment of the signs 6j = ± in the expression of the tangent vector to 
each link: ejXj + CjXj. Due to the folding rule (|1.2|), each Xj must appear with the two 
signs, therefore only three signs have to be specified for each triangle, leading to 2 3 = 8 
possibilities. 

We recover in this way the 48 = 6 x 8 face configurations mentioned above. This 
presentation has the advantage of decoupling the colouring step (i) from the assignment 
of signs (ii), therefore displaying the degeneracy of the 3d folding when compared to 3- 
colouring. 

For a given 3-colouring of the links of the triangular lattice, let us now evaluate 
the number of choices of signs we can make. This is best done on the dual hexagonal 
lattice, whose vertices are the centers of the faces of the triangular lattice, and whose 
links (represented by dotted lines in Fig. 12) cross those of the triangular lattice. The 
tangent vectors are now link variables of the hexagonal lattice. A 3-colouring of the links 
of the triangular lattice is simply a colouring of the links of the hexagonal lattice, with the 
constraint that the three links meeting at each vertex have distinct colours. Now a choice 
of sign, say e\ for the coefficient of xi has to be made on all B and R links. Moreover, this 
sign propagates from a given B (resp. R) link to the two neighbouring R (resp. B) links, 
in order for ( |1.2| ) to be satisfied at each vertex. 

More precisely, as shown in Fig. 13, the sequences of £?, R, B, R, ... links form closed loops 
on the hexagonal lattice, and the sign of the xi component must alternate along each such 
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I 

W B W 

W^B^R^W R^B^R W 

R B W B 

W B ^ R ^ W R ^ B ^ R ^ W 

W B W 

I 

Fig. 13: A sample 3-colouring of the hexagonal lattice: the BRBR... se- 
quences of links form dense loops, represented by thick lines. The coefficient 
of xi in the corresponding tangent vectors can be independently fixed on each 
loop to be either +1 on the B links and —1 on the R links, or —1 on the B 
links and +1 on the R links. 

loop. For each BR loop, there are two choices of signs of the Xi components: e\ = + (resp. 

— ) on all B (resp. R) links of the loop, or conversely. The same holds independently for 

the choices of signs of the X2 components along the BW loops and of the X3 components 

along the WR loops. 

So for a given colouring of the links of the hexagonal lattice, we are left with 2 choices 
of signs for each of the mi BR loops, m<i BW loops and WR loops. The partition 
function for 3d folding reads then 

Z 3d = 2 mi+m2+m3 (5.1) 

3— colourings 

where the sum extends over all the 3-colourings of the links of the hexagonal lattice and mi 
(resp. m2, friz) denote the number of BR (resp. BW, WR) loops for each colouring. The 
expression ( |5.1| ) identifies the 3d folding problem of the triangular lattice with a dressed 
3-colouring problem of the hexagonal lattice, obtained by attaching a Z2 variable e to 
each loop of any alternating two colours. Notice finally that the bi-coloured loops above 
form three dense coverings of the hexagonal lattice, in the sense that each vertex of the 
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hexagonal lattice belongs exactly to one BR, one BW and one WR loop. The loops of a 
given type are non-intersecting. 

As an exercise let us re-derive the number 96 of vertex configurations in this dressed 
3-colouring framework. This number is simply the partition function of a single hexagon 
H with external legs (dual to an elementary hexagon of the triangular lattice) 

Z H = 2 mi+m2+m3 . (5.2) 

3— colourings of H 

Two situations may occur: 

(i) the colouring contains a bi-coloured central loop, in which case one of the m's is equal 
to 1 and the two others are equal to 3 (open loops), leading to a weig ht 2 7 . There are 6 
such colourings. 

(ii) all the bi-coloured loops are open. Each external leg belongs to two loops and each 
loop contains two external legs; the total number of loops is thus equal to the number of 
external legs 6: this leads to a weig ht 2 6 . There are 6x11 — 6 such colourings (since the 
total number of 3-colourings of H is 6 x 11). 

We finally get 

Z H = 6 x 2 7 + 60 x 2 6 = 4608 (5.3) 



in agreement with (|3.4|) . 



5.2. Simple bounds on the entropy 

The simplest lower bound on q^a = limjv A -Kx> z\d A i ls obtained by minoration of the 



sum (|5.1|) by picking a particular 3-colouring of the hexagonal lattice and evaluating its 
contribution. 
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f^m w r b t^m 

L l 

W^B^R W B R W B 

B R W B R W B 
B 



WB RWB RWB 
K B W R 



Fig. 14: The antiferromagnetic groundstate which maximises mi, W2 and 
W3 simultaneously. The Si? (resp. BW, VFi?) hexagonal loops are indicated 
by discs of different colours. 

We take the particular configuration depicted in Fig. 14, which maximises mi, m2 and 
ms simultaneously. This configuration corresponds to a regular antiferromagnetic ordering 
of the a variables. The B and R links are arranged in a dense regular set of hexagonal 
loops, as well as the B and W links and the W and R links. Since the smallest loop on 
the hexagonal lattice has length 6, mi, m,2 and ms are clearly maximal. Noticing that 
the total number Nh of hexagons in the hexagonal lattice is equal to the total number 
N v = of vertices of the dual triangular lattice, this configuration has 

Nh N v N a 
mi = m 2 = m 3 = — = — = — . (5.4) 

Hence Z 3d > 2 Na / 2 , so that 

q 3 d > V2 = 1.414... (5.5) 

In the 96-vertex language, the contribution of the antiferromagnetic groundstate of Fig. 14 
corresponds to a restriction of the model to the subset of vertices with alternating a 
variables. This subset is characterized by the presence of a thick or a thin solid line on 
each internal link (see Fig. 10). There are 1 + 6 + 3 + 6 = 16 such vertices. In this respect, 
\pl is simply the exact partition function per triangle of this restricted 16-vertex model, 
as can be checked directly. 
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We can improve this bound by incorporating more configurations in our counting, i.e. 
excitations of the antiferromagnetic groundstate of Fig. 14. This is done in the Appendix, 
leading to the improved lower bound 

q 3d > 1.429... (5.6) 

Using now the fact that, for an arbitrary configuration of 3-co louring of the hexagonal 
lattice, the numbers m±, m2 and m 3 are always smaller or equal to N&/6, we also get a 
simple upper bound 

Z 3d < 2^/ 2 1 = ^ Z ™ ( 5 - 7 ) 

3— colourings 

since the 2d folding problem and that of 3-colouring of the hexagonal lattice are equivalent 
0. Thus, we get an upper bound on q 3d in terms of q 2d ( P--3Q 

q 3d < V2q 2d = 1.709... (5.8) 

5.3. Improved bounds from 2d folding in a field 

The partition function of the 2d folding problem also has an expansion in terms of 
dense loops on the hexagonal lattice. As mentioned above, the 2d folding problem is 
equivalent to that of 3-colouring of the links of the hexagonal lattice, dual to the original 
triangular lattice. Instead of considering the 3-colourings of the hexagonal lattice links, 
however, we can concentrate on say the VF-coloured links only. Consider a particular 3- 
colouring of the links. We have seen above how the paths of BRBR... links form dense, 
non-intersecting closed loops on the hexagonal lattice. Exchanging the B and R links 
along any of these loops independently leads to equally admissible 3-colourings of the 
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links. Therefore, for a given admissible configuration of W links, one is left with 2 possible 
independent choices of colourings per BR loop. The number mi of such loops is fixed by 
the position of the W links only. We can write 

z ** = E 1 = E 27711 ■ ( 5 - 9 ) 

3 — colourings W link config. 

This is the 2d version of the 3d folding dense loop expression (|5TT|). 

The solution by Baxter || of the 3-colouring problem of the links of the hexagonal 
lattice includes the introduction of an extra parameter in the weighting of configurations. 
This parameter may be interpreted as a staggered magnetic field h s t in the following way. 
In the language of the face spin (a) variable (a = +1 or —1 on a face indicates whether 
the colours of adjacent links around the face increase (B — > W ^ R — > B) or decrease 
(B — > R — > W — > B) counterclockwise), this weight reads 

e ahst per triangle facing up, A , 

(5.10) 

e -o- st p er triangle facing down, V . 
This gives a maximal weight to the completely folded (antiferromagnetic) groundstate of 
the model (see Fig. 14). Baxter's result is the exact partition function per triangle 



'1 _ 2 12(l-3nV 
2,12(2-3™)) (1 _ z 12(-3n)) 



q2d(h st ) = z H y/1 \Z„ ' ^ 5 ^ ( 5 - n ) 

n=l V I 1 



with z = e hst . The abovementioned result ( p..3| ) for the 2d folding entropy is then recovered 
in the limit h s t — > 0, i.e. z — > 1. 

Let us now translate the result (|5.11|) into the language of dense BR loops (|5J^) on 



the hexagonal lattice. Consider a particular admissible configuration of W links (a term 
in the sum in Q5.9Q ). The extra weights ( |5.10| ) can be translated into a dressing of the 
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lt\ A ^ - 5 



Fig. 15: The four basic possibilities of local orientation for the succession of 
a B and a R oriented link on the hexagonal lattice (up to global rotation by 
±27r/3). The corresponding state of the visited (dual) triangle is represented 
above each two-link state. The weight in a staggered magnetic field is e hst in 
the first and third case, and e~ hat in the second and fourth ones. The former 
correspond to left turns, whereas the latter correspond to right turns along 
the BR loop. Exchanging the B and R links would simply reverse the sign 
of the face spin a as well as the orientation of the links. 

corresponding BR loops. More precisely, a succession of two B and R links along a BR 
loop may be in one of the 4 states (up to a rotation by ±27r/3) depicted in Fig. 15. This 
provides us with a natural orientation of the loops, according to Fig. 15. Each link joins the 
centers of two adjacent triangles, one pointing up and the other pointing down. We orient 
the B links toward the upward pointing triangle, and the R links toward the downward 
pointing triangle. Since triangles pointing up and down alternate along a loop, as well 
as B and R links, this rule defines a consistent orientation for each BR loop. Reversing 
the orientation of a loop amounts to exchanging its B and R links, with fixed W links. 
By inspection of the 4 cases of Fig. 15, we see that the weights ( p.!0| ) translate into local 
weights 

e hst per left turn, 

(5.12) 

e~ hst per right turn 

of each oriented BR loop. The previous degeneracy of 2 per loop now amounts to summing 
over the two orientations of each loop, resulting in a total weight 

e 6hst + e~ 6hst = 2cosh6/i st (5.13) 
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per loop, as a closed loop has its number of left turns minus its number of right turns 

equal to ±6, and the opposite for the opposite orientation. The partition function for the 

2d folding problem in a staggered magnetic field reads therefore 

Z 2d {h st ) = (2cosh6/i st ) mi 

W link config. (5-14) 

oc q 2d (h st ) NA 



where the sum is as in ( |5.9| ) and q2d(h s t) is defined in ( |5.11| ). To get a more symmetric 



form, we can sum as well over the 3-colourings of the links of the hexagonal lattice, thus 

removing a factor 2 mi in the above expression, hence 

Z 2d (h st ) = (cosh6/i st ) mi x2 mi 



W link config. 

= E (coshG^rE 1 (5.15) 

W link config. B,R 

£ (cosh6/i st ) mi . 

3— colourings 

In the last expression, we could replace mi by ra 2 or 1713 without distinction. 

With this function at hand, we are now ready to give more bounds on q^a- Starting 
from a given 3-colouring configuration of the dense loop expression ( |5.1[ ), let us examine 
the numbers m±, m 2 and in more detail. A given BR loop is called direct (resp. 
indirect) iff its total weight (|5.12| ) is e Ghat (resp. e~ 6hst ). We can therefore decompose 



the number mi of BR loops into the numbers d± and i\ of direct and indirect BR loops 
respectively. Analogously, we can define the numbers d 2 , ds, i 2 , 13 of direct and indirect 

BW and WR loops. We have the following relations 

m j = dj + ij , j = 1, 2, 3 

(5.16) 

d\ — i\ = d 2 — % 2 = d 3 - is . 
To prove the last identity, we note that for a fixed 3-colouring, the total weight ( |5.10| ) 

can be evaluated in terms of any of the three systems of dense loops (BR, BW or WR) 

32 



independently, with the same result 

e ® h st(J2 A = e Gh st (di-ii) _ e <5h st (d 2 — h) _ g6/i s t(d 3 -i3) (5-17) 

A first lower bound on can be obtained as follows. In the dense loop expression 
( p . 1| ) of the partition function Z^-, we restrict the sum to the configurations where all the 
BR loops are direct. For each admissible colouring of the W links, there is exactly one 
such configuration of B and R links (it corresponds to a particular choice of orientation 
on each loop). By summing over these states, we get a minoration Z^ of Z 3 d. All these 
states have by definition i\ = 0, hence mi = d x . Using the identities ( |5.16| ), we find that 



(5.18) 



m 2 = d 2 + i 2 = {d 2 - h) + 2^2 = d x + 2i 2 = mi + 2i 2 > mi , 
^3 = d 3 + is = (c?3 - ^3) + 2z3 = di + 2i 3 = mi + 2i s > mi . 
The partition function can therefore be bounded from below by 

Zj* r R = 2 mi+rn2+ma > Y 23mi • ( 5 - 19 ) 

W link config. W Jink config. 

all BR loops direct 



The latter is readily identified as the partition function (|5.14|) of 2d folding in a staggered 
magnetic field h* t , such that 

2cosh6/i* t = 2 3 (5.20) 



i.e. 



h* st = (4 + Vl5) . (5.21) 



Finally we get the bound 



q 3 d > q d B r R > q2d(Kt) = 1-421... (5.22) 
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If for a given W link configuration, we now take into account all the orientations of 



the BR loops, eq. (|5.18| ) no longer holds, but we still have the inequalities 

m 2 = (rf 2 -«2) + 2z 2 = (di—ii) + 2i 2 = mi + 2i 2 - 2i\ > m\—2i\ 

m 3 = (d 3 - i 3 ) + 2i 3 = (di — it) + 2i 3 = rri\ J r2i 3 -2i\ > m\—2i\ 
so that we can write 



(5.23) 



z 3d > Yl 2 3mi - 411 = 23mi Y 2 ~ 411 ( 5 - 24 ) 

3 — colourings W link config. B,R links 

where we simply separate the sum over the admissible W links from that over the B and 
R links along BR loops. Since the mi BR loops can be chosen to be direct or indirect 
independently, the latter sum factorises and contributes for 

mi / • \ i i 

E I 2^ = (! + ^) mi = (17/16)-, (5.25) 
and hence we finally get 

Z 3d > Y (l 7 / 2 )" 1 • ( 5 - 26 ) 

W links 

According to eq.(p.l4|), the latter partition sum is nothing but that of 2d folding in stag- 



gered magnetic field h*i , such that 



17 

2cosh6/i* t * = — (5.27) 



i.e. 



6 

This gives the following lower bound on q 3d 



= 1 17W273 _ (Mg) 



qsd > q2d(K* t ) = 1.4351799... (5.29) 
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Comparing this bound with the numerical estimate (|4.8|) , we suspect that the exact 
result for q^d differs from this value by no more than 1 percent. 

We can also obtain a majoration of Z%d by use of the Holder inequality on averages, 
namely 

(AxBxC) < (A Q )«(S /3 )^(C 7 )^ (5.30) 

where 1/a + 1/P + l/'j = 1. The majoration reads 

Z 3d = 2mi x 2 m2 x 2 ms 

3— colourings 

< ( j2 2 ami I x 

\ 3— colourings 



(5.31) 



x I > 2 /3ma I x 



£ 21 

3— colourings 



T 



X 



2 7ma 

3— colourings 

The 3 terms on the rhs of (|5.31| ) are readily identified as powers of 2d folding partition 



sums ( pM4| ) with respective staggered magnetic fields 



h st (x) = hn (2 X + V2 2 *-1) x = a,/3, 7 . (5.32) 

D 

The lowest upper bound provided by (|5.30|) corresponds in fact to a = (3 = 7 = 3, which 
leads to 

Q3d < g 2d (M3)) = 1-589469... (5.33) 

6. Discussion 

6.1. Estimating q^d 

In our strategy for getting lower bounds on ^3^, we first considered the contribution 
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of a particular state, the groundstate of Fig. 14. In this state, the W links form a regular 
hexagonal pattern, and the mi = BR loops are all direct. This leads to the first 

estimate v^2 of eq.( |5.5|) . Keeping the W links fixed, elementary excitations are obtained 
by exchanging the B and R links along some loops, thus reversing their orientation from 
direct to indirect. As shown in the Appendix, these excitations have a fugacity 1/2 4 
along with some interactions. Ignoring the interactions leads to the improved lower bound 
V / 2(l + 1/2 4 ) 1 / 6 (A. 3). Taking these interactions into account leads to the lower bound 
Qbr of (A. 5), expressed in terms of the entropy of a loop gas on the hexagonal lattice, 
reproducing the high temperature expansion of the 0(n = 4) model at coupling K = 1/2 
(A.8). 

In parallel, we also considered the sum over all admissible configurations of W links. 
With all BR loops direct, we get the first estimate q2d(h* t ) of eq.( |5.22|) , with 2 cosh6/i* t = 



\f 7 t ' . For each W link configuration, elementary excitations also correspond to reversing 
the orientation of some BR loops from direct to indirect. There also, we were led to 
assign a weight 1/2 4 per indirect loop (see (|5.24|) ), and we obtained the bound q2d{h*l) of 



eq.( |5.29|) with 2 cosh6/i*^ = (\/2(l + 1/2 4 ) 1 / 6 ) 6 . It is therefore tempting to conjecture the 
lower bound 

qsd > thdiK?) (6-1) 



with 



2cosh6/C** = {qbr) 6 (6.2) 



to account for interactions between the indirect loops, with qsR given by (A.8). This 
lower bound could in fact be the exact value of q^a since it would incorporate all the 
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effects of excitations, qsn is estimated through a Mayer expansion in eq.(A.4). With this 
estimation, we get 

q2d(hT) ~ 1-4356... (6.3) 

6.2. d-dimensional generalisation 

The 3d octahedral folding problem has a natural generalisation to R d , by considering 
folding on a generalised <i-dimensional FCC lattice. Equivalently, this is done by restricting 
the (i-dimensional tangent vectors to be the edges of oriented polytope of H d generalising 
the octahedron, whose 2 x d vertices have positions ±ei/v^2, where e$, % = 1, d is the 
canonical basis of R d . The 2d(d — 1) unit tangent vectors read e^x^ + e^-Xj (xj = e^/y^ 
as before). They form 4d(d — l)(d — 2)/3 triplets with vanishing sum, corresponding to 
the faces of the polytope. This gives 3! x Ad{d - l)(d - 2)/3 = 8d(d — l)(d — 2) possible 
environments for a given triangle. Any such environment still takes the form displayed 
in Fig. 12, where now l<i^j^k<d. As in the 3d case, let us consider the tangent 
vectors as link variables on the dual hexagonal lattice. The tangent vectors with non-zero 
Xj component form rrii loops along which the sign q alternates, for i = l,...,d These 
loops form a dense covering of the hexagonal lattice in the following way: any vertex of 
the hexagonal lattice belongs to exactly three loops of different type, and any link of the 
hexagonal lattice belongs to exactly two loops of different type. As in the 3d case, the two 
choices of signs per loop lead to a partition function 

Z d = 2 mi+m2+ - +md . (6.4) 

dense loops 

To compute the number of possible vertex environments, we evaluate the partition 
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function of a single hexagon H with external legs (as in (|5.2|) ) 



miH \-m d 



dense loops 
on H 



(6.5) 



2'r d + 2 b s d 



where r d (resp. Sd) denotes the number of configurations with (resp. without) an internal 
closed loop: the first term corresponds to 7 loops (1 closed internal loop and 6 open loops 
entering and exiting the hexagon through the external legs), whereas the second only has 
a total number of 6 loops. Let us first compute the number of configurations with a 
closed internal loop. We fix the type of the central loop, say to xi (among the d possible 
choices). Around each vertex of the hexagon, as the xi loop occupies the two inner (left 
and right) links, the only possibility is that a loop of type x^ occupies the left and external 
links, whereas a loop of type Xj occupies the right and external links, with 2 < i ^ j < d. 
This suggests the introduction of a transfer matrix of size (d— 1) x (d— 1) mapping 
an internal link to the subsequent one, with entries 



where the prefactor d accounts for the d possible choices of the type of the internal loop. 
Next let us evaluate the total number + Sd of loop configurations on H. Here again, we 
need to construct a transfer matrix from an internal link to the subsequent one. The state 
of a link is specified by the two types i < j of loops to which it belongs, namely d(d — l)/2 
choices. The desired transfer matrix has therefore the size d<yd 2 1 ~ ) x d Z > an< ^ its entries 




(1-<M 2<i,j<d. 



(6.6) 



We find 




(6.7) 
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read 

Mij-ki = 6 ik (l - 8ji) + 8ji(l - 8 ik ) + 6u + 8 jk (6.8) 

for 1 < i < j < d and 1 < k < I < d. Note that this matrix is slightly different from that 
used for the 3d folding problem (|3~3|). After some algebra, we find 



r d + s d = TrM 6 = d(d - l)(d - 2)(d 4 + 42d 3 - 380d 2 + 1096d - 1072) . (6.9) 



Using (|6.7| ) and (|6.5|), we finally get 

Z djH = 64d(d-l)(d-2)(2d 4 + 33cf 3 -349d 2 + 1047d-1041) . (6.10) 

Factoring out the 8d(d — l)(d — 2) configurations of a given triangle, the (i-dimensional 
folding problem is therefore equivalent to a V^-vertex model, with 

V d = 8(2d 4 + 33d 3 - 349d 2 + 1047d - 1041) . (6.11) 

We recover V3 = 96 for the 3d folding problem, whereas V4 = 1496, V5 = 6752, Ve = 
19176... 



7. Conclusion 

In this paper, we have defined the 3d octahedral folding i.e. the folding of the triangu- 
lar lattice on the 3d FCC lattice. This model was formulated as a 96-vertex model, with 
two face spin variables z and a subject to the two folding rules (|3.9|) ( |3.13|) . Equivalently, 
the partition function of this model was reexpressed as that of a dressed 3-colouring prob- 
lem, involving a dense covering of the hexagonal lattice by bi-coloured loops. With these 
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two formulations at hand, we were able to estimate the folding entropy ssd = hi q^d, both 
numerically by use of a transfer matrix and analytically by deriving various exact bounds. 

Beyond mere counting of folding states, it would be interesting to obtain the complete 
phase diagram of this system, including both bending rigidity and magnetic field as was 
performed in |J for the 2d case. In particular, it would be desirable to know the precise 
status of the crumpling transition in this framework, including its order (first or second). 
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Appendix. Lower bounds on the entropy from local excitations 




Fig. 16: An elementary BR excitation of the antiferromagnetic groundstate 
of Fig. 14. The B and R links are exchanged on the shaded hexagon. The 
three neighbouring BW hexagons are glued to form one loop (thick solid line 
on the first figure), whereas similarly the three neighbouring WR hexagons 
are glued to form one loop (thick dashed lines on the second figure). 

Starting from the fundamental colouring state of Fig. 14, we consider the following local 
excitation: pick a particular BR hexagon, and exchange its B and R links. The effect of a 
BR excitation on the neighbouring BW and WR loops is illustrated in Fig. 16. In this new 
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colouring configuration, mi is clearly unchanged, but the three BW neighbouring hexagons 
have been glued into just one loop (thick solid line on Fig. 16), hence is decreased by 2, 
and the same applies to the three WR neighbouring hexagons (thick dashed line in Fig. 16), 
hence ms is also decreased by 2. Therefore the total weight of the excited configuration is 



2 2 . 



(A.l) 



Now there may be any number of BR local excitations in the system, but their interactions 
will result in a modification of their weight. The number of loops can only be increased by 
such interactions however, as will become clear below. From now on, we factor out the term 
\/2 A in the partition function. The above leads to a lower bound for Z^d by considering 
a non-interacting gas of BR excitations with fugacity 1/2 4 , taking place among the Nh/3 
hexagons of BR type, with partition function 

~(o) _ ^ !^ 



, BR - 2— x(l + -)^ 



17 

T 



N A 



Hence the improved lower bound on q 3 d 



Q3d 



> 



17 \ « 
T 



1.4285. 



(A.2) 



(A.3) 




Fig. 17: Two neighbouring BR excitations. The BR excitations take place 
on the shaded hexagons. The first figure shows how five former BW hexagons 
have been glued into a single loop (thick solid line) leading to — > ni2 — 4. 
The second figure shows how five former WR hexagons have been glued into 
a single loop (thick dashed line) leading to 1713 — > 7713 — 4. 
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The BR excitations take place on the BR hexagons of the fundamental state of Fig. 14, 
i.e. on the vertices of a triangular lattice of larger size, which we call the excitation lattice 
from now on. Excitations which are not neighbours on this lattice do not interact, namely 
their total weight factorises into the product of their individual weights. In fact, there 
is no two-body interaction between them: when two such excitations are neighbours as 
illustrated in Fig. 17, i.e. only connected by a W link, the total number of loops is reduced 
by 8 (five neighbouring BW hexagons are glued into one loop, hence 1112 — > ni2 — 4, and 
five neighbouring WR hexagons are glued into one loop, hence — > rris — 4); hence a 
relative weight 1/2 8 = 1/2 4 x 1/2 4 . 




Fig. 18: Three neighbouring BR excitations. The picture on the left shows 
how the 6 former neighbouring BW hexagons have been glued to form two 
loops (thick solid lines) leading to mi — > — 4. The picture on the right 
shows how the 7 neighbouring WR hexagons have been glued to form a single 
loop (thick dashed line) leading to — > 7713 — 6. The total number of loops is 
therefore decreased by 10 = 12 — 2, showing that the three-body interaction 
weight is 2 2 . Note the asymmetry between BW and WR loops. 

There is, however, a non-trivial three-body interaction, when three BR excitations 
take place on the three BR neighbours of a BW or a WR hexagon, as depicted in Fig. 18. 
Indeed, the three excitations re-create a central loop with B and R exchanged, thus 
increasing the total number of loops. A more precise counting (see Fig. 18) shows in this case 
that the total number of loops is increased by 2 (only 10 loops have disappeared, instead 
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of 12 = 3 x 4 for three independent excitations), resulting in a three-body interaction 
weight 2 2 . The partition function Zbr incorporating any BR excitation can be evaluated 
order by order in the number of BR excitations. At this stage, the lower bound (|A.3| ) 
can therefore still be improved by performing a Mayer expansion of the free energy per 
hexagon, in increasing order of the number of BR excitations. Up to order 5 in the number 
of BR excitations, we find the partition function per hexagon 

1 3 3 69 
(q BR ) 6 = 8(H^ + ^ + ^ + ^ + -)> 8 x 1.0644... (A.4) 

hence the improved lower bound 

qsd > Qbr > 1-429... (A.5) 

Let us finally show that the contribution Zbr to the partition function of the most 
general combination of BR excitations of the above groundstate is directly expressible as 
the high temperature expansion of a particular 0(n) model on the hexagonal lattice. 
Indeed, any combination of excitations can be decomposed into clusters which do not 
interact with each other. A typical cluster of BR excitations is shown in Fig. 19: the BR 
hexagons are represented by thin dotted lines; the excited BR hexagons (i.e. with B and R 
links exchanged) forming the cluster are shaded. Note that the W links stay in the initial 
state of Fig. 14, as they are not affected by BR excitations. The clusters are indeed clusters 
of neighbouring vertices on the excitation lattice. The dual white links form an hexagonal 
sublattice of the original triangular lattice, which is also the dual of the excitation lattice. 
We define the boundary of a given cluster as its envelope of W links on this lattice (the 
boundary does not include internal W links which separate neighbouring BR excitations in 
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Fig. 19: A typical cluster of BR excitations. The BR hexagons are repre- 
sented by thin dotted lines: their centers sit at the vertices of a triangular 
lattice. The shaded hexagons are excited, namely their B and R links are ex- 
changed. The two-component thick solid line is the boundary of the cluster, 
drawn on the links of the original triangular lattice (thin dotted lines) . The 3 
loops in thick dashed lines result from the gluing of the BW hexagons touching 
the boundary of the cluster. This reduces number of BW loops m 2 — > m 2 — 10. 
The WR loops have not been represented for simplicity. The reader will con- 
vince himself that their number is reduced by the cluster to 777.3 — > wi3 — 12. 
On the other hand, the two connected components of the boundary have re- 
spective lengths 6 and 20, satisfying (6 — 2) + (20 — 2) = 10 + 12, as expected. 

the cluster). The boundary of the cluster is represented in Fig. 19 by thick solid lines: note 

that the boundary links (which are W links of the original triangular lattice) cut white 

links of former BW or WR hexagons, which are glued together by the BR excitation. A 

given cluster C receives the relative weight 

n 22-4 ( a - 6 ) 

boundary connected 
components c 

where the product extends over all the connected components c of the boundary of the 
cluster (there are two such components for the boundary of the cluster of Fig. 19). The 
number l c denotes the total length of the connected component c, namely the total number 
of W links forming c (the two corresponding lengths are 6 and 20 for the example of Fig. 19) . 
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This result can be easily understood qualitatively. In the bulk inside a cluster, the 
hexagonal loops of BW and WR type are simply exchanged. Therefore, the changes in 
the numbers of loops can only take place at the boundary of the cluster, where some BW 
(resp. WR) hexagons are glued together to form larger loops, thus decreasing the values 

of 777.2 (resp. 777,3 ). 

Eq. (|A.6|) can be proved by recursion, by studying the effect of adding an extra BR 



excitation to a given cluster C. Of course ( |A.6|) holds for a single isolated BR excitation 



in which case l c = 6, and we recover the individual weight 1/2 4 ( |A.1| ). 

Notice now that the weight (|A.6|) only depends on the boundaries of the clusters, 



which form a set of loops drawn on the white links of the groundstate of Fig. 14. This leads 
to a reexpression of the partition function Zbr for all possible BR excitations as that of 
a loop gas on this (original W links) hexagonal lattice 



'BR 



2 ^ E n 22 "' ( a - ? ) 

loops on the loops 
hexag. latt. 

where I denotes the length of each loop. The sum in eq. (|A.7|) is nothing but the high 



temperature expansion of the 0(n) model of |12] on the hexagonal lattice, with n = 4 and 
K = 1/2. Hence 

Qbr = V2(q oin =A)(K = l/2)y 6 (A.8) 

where qo( n ){K) denotes the thermodynamic partition function per hexagon of the 0(n) 
model. Unfortunately the partition function for the 0{n = 4) model is not known exactly. 
The model can be mapped |H| onto a 6-vertex model on the Kagome lattice, whose vertices 



sit at the center of the links of the hexagonal lattice. The latter however is a 6-vertex 
model in a magnetic field, for which no exact solution is available. Nevertheless, the Mayer 
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expansion above ( |A.4|) is expected to converge rapidly, and the value found for the lower 
bound should be accurate to the first three digits (the order 5 term in the expansion is 
only ~ 10" 4 ). 

Still qsR as given by (|A.8| ) is only a lower bound on q^d- 
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